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Abstract 

Some simple (namely, single-channel) correlation functions involving an arbitrary number 
of fields are computed by means of a direct application of the residue calculus, through 
partial fraction expansions. Examples are presented in minimal models and parafermionic 
conformal theories. A generic factorized expression is deduced for the corresponding single- 
channel structure constants. 
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1. Introduction 

1.1. Correlators, OPEs and Ward identities 

The plain method for evaluating Appoint correlation functions, given the OPEs (op- 
erator product expansions) of the fields whose correlations are to be computed, amounts 
to substitute the OPE of two fields in order to reduce the correlator to (N — l)-point 
functions and iterate this procedure until the result becomes expressed in terms of a three- 
point function. The later being known exactly up to structure constants, the correlator 
is then expressed in terms of these constants. However, the reduction in the number of 
points has been traded for a new complication: we then have to sum up the infinite series 
(i.e., conformal blocks) associated to each OPE. In principle, these can be summed exactly 
only in simple cases. A sample computation is presented in Appendix A. 

For the mere formulation of the model, the correlators that are particularly important 
are those involving the symmetry generators, namely the generators of the extended con- 
formal algebra. Their relevance lie in that the internal coherence of the extended algebra 
boils down to a precise statement concerning the correlators of the symmetry generators: 
all their four-point functions must be associative. The associativity requirement is the con- 
dition that a correlation function can be calculated in many different ways, in particular, 
by evaluating the OPEs in different orders, without affecting the result. (This is the way 
the structure constants are calculated.) But testing associativity calls for the exact form 
of these correlation functions. 

However, in this particular instance, the problem appears to be tractable. Indeed, 
the evaluation of correlation functions involving extended conformal-algebra generators 
is usually rather simple in that only the singular terms have to be considered. Take 
for instance a correlation function involving the energy-momentum tensor T(z) and some 
primary fields FT/Li ^j( z j)- This is certainly relevant to the question of studying the 
associativity of a given conformal algebra because the algebra generators other that T 
have to be Virasoro primary fields. It is a very basic fact that in eliminating T(z) through 
OPE, one simply needs to take into account the singular terms in the product of T(z) with 
all the other fields in the correlator: 



In other words, in this special case we do not have to keep track of an infinite series. This 
is a consequence of the conformal Ward identities - cf. [1, 2, 3]. 
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1.2. Correlators involving T as meromorphic functions 



But if we think about this result from the point of view of OPEs, it looks rather 
surprising that by considering only the singular terms of the OPE of T with the other 
fields of the correlator (and not the complete infinite OPE series) we can compute the 
correlation function exactly. This has a natural complex-analysis explanation: the result 
simply corresponds to the partial fraction expansion of the meromorphic function repre- 
senting the correlation, viewed as a function of one of its field. For instance, the correlator 
in (1.1), considered as a function of z, is a meromorphic function with double poles at the 
various Zj (since the 0j(zj)'s are supposed to be primary) with coefficients fixed by the 
OPE T(z)<pj(zj). In that case, the partial fraction expansion is complete in that there is 
no additional analytic piece. Indeed, this meromorphic function vanishes at infinity since 
T(z) ~ z~ 4 as z^oo. Such a function is simply given by the sum of the principal parts at 
the various poles. 

Let us make the above statements more explicit. Recall that a meromorphic function 
with vanishing analytic part can be written in partial fraction as 

JV «j (j) 

j=lr=1 y z z j> 

where nj is the order of the pole at Zj and the coefficients a,r are given by 

^-'-.^H^-^W (L3) 

When F(z) is a correlation function, the various coefficients are expressed in terms of 
lower-order correlation functions. On the other hand, if the meromrophic function F(z) 
does not vanish at infinity, an analytic function needs to be added to this sum of principal 
parts. Suppose that F(z) behaves rather like F(z) ~ z p as z^oo, with p integer. This 
signals the presence of a pole of order p at infinity so that the principal part at infinity (in 
which we include a constant term) has to be taken into account. The expression of F(z) 
becomes then 

j=l r=l ^ i' k=l 

The last sum is the analytic part of F(z). 



It is then completely obvious that in the case (1.1), the principal part at Zj is nothing 
but the sum of the singular terms in the OPE T(z)(j)j(zj). Indeed, viewing the correlator 
as a function of z, we have: 



j — 1 i — 1 \ 



with 



N 
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= a i (IJ^(^)> 
J = l 

(the first term of the second line drops because it is independent of z); this gives the rhs 
of (1.1). 

The above derivation of (1.1) makes clear the rather auxiliary aspect of the primary 
nature of the fields inside the correlator. For instance, one could introduce quasi-primary 



fields, e.g., T(z) itself. Regarding the following correlator as a function of £1 leads to [2]: 

M N 

<n T ^)n<^-)> 

i=i j=i 

M M N 

= E 77^1(11^)11^^)) 

i=2 ^ ^ fc=2 j = i 

fe#l (1.8) 

A manifestation of the associativity property is that the same expression for this correlator 
follows by considering it as a meromorphic function of another variable Q. 

Even 'less- primary' fields are treated by the same method: only the singular terms of 
the OPE with T(z) contribute. 

1.3. Generalizing the class of correlators computable from partial fraction expansions 

What is central in the previous computations of correlators involving T? At first, 
there is the fact that the OPE of T with another field has only one channel, i.e., all those 
terms that appear in the OPE T(z)<p(w) belong to the conformal family of <p{w). This 
ensures that the different powers of z — Zj (for a given j) all differ by integers. A second 
important point is that T is local with respect to any other field in the theory, which 
implies that, not only the different powers of z — Zj differ by integers, but they are all 
integers themselves. In other words, the OPE T(z)(fi(w) is a genuine Laurent series. 

Therefore, this computation method can be used for calculating the different four-point 
functions for the extended-algebra generators when these generators X all have integer 
dimensions. Note that in this case, the one-channel constraint is superfluous because all the 
powers of z — w in the OPE yi(z)yj(w) G X]/ C [^fc(' u; )] have automatically integer dimension. 
A well-known example of this type is the WZW model whose extended symmetry is an 
affine Lie algebra, with generators J a satisfying the OPE J a x J b = I + J c . 

But if this is to be used for the analysis of more general extended algebras, for instance 
parafermionic models with fields having fractional dimension, one has to be able to tackle 
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situations where there are branching points, that is, when the powers of z — w are fractional 
in the OPE yi(z)yj(w). The cure is quite simple: we just have to modify the correlation 
function, which is viewed as a function of the first variable z\ - the position of a given 
symmetry generator, say y{z\) -, by multiplying it by appropriate powers of z\ — Zj (zj 
being the position of another field inside of the correlator) in order to transform it into a 
meromorphic function of z\ [4]. In other words, to analyze the correlation (y(zi) ■ ■ •), we 
consider the intermediate function 



for those values of di appropriate to make F(zi) meromorphic. The OPEs fix then the 
position of the poles of this function together with their residues. This determines F(zi) 
and, thereby, the correlation function under consideration. Notice however that in the 
present case, the principal part of F(z±) at Zj is not given solely by the sum of the singular 
terms in z\ — Zj\ for the subleading terms in the principal series, there are derivatives in the 
expression for the constants a^_ k in (1.3) and these derivatives do not select exclusively 
the corresponding singular terms due to the presence of the prefactors Yli>2( Zli ) di on 
which the derivatives also act. 

As already pointed out, the presence or absence of a regular (analytic) part in the 
expression of F{z\) is fixed by the behavior of the field at position z±, as zi^oo, e.g., 
~ Zi y , together with the prefactor composed of the different fractional powers of 
z\ — Zj that have been introduced. A modification of the prefactor by an integer power of 
z\ j obviously affects the nature of the regular part of F(zi) but not the final expression of 
the correlator. 

By bringing out the conditions underlying the applicability of the residue method 
to the computation of correlators involving T, we have identified at once two criteria: 
(1) the single-channel requirement and (2) the locality condition. The simple trick just 
described for transforming a function with branching singularities into a meromorphic 
function provides thus a way of bypassing the apparent second limitation. 

In relation with the single-channel requirement, we have already presented a situation 
in which it can be relaxed, namely for correlators of symmetry generators all having integer 
dimension. Phrased in more general terms, the one-channel condition is not mandatory 




(1.9) 
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when the OPEs of the fields inside the correlators close in a set of fields which all have 
conformal dimension that differs from each other by integers. But this is a rather special 
instance. In more general circumstances, the single-channel requirement cannot be avoided. 
Notice however that this is not as restricting as it looks at first sight. What is really needed 
is a not exactly a genuine single-channel OPE but, rather, a correlation function that selects 
a single channel in each intermediate OPE. 

In the above considerations there is in addition an implicit third limitation, which is 
that the field in terms of which the partial fraction is formulated has to be a symmetry 
generator. But the complex analysis is blind to the subtle conformal nature of the field 
evaluated at z\ . It is clear that it can be any field, as long as (generically) the correlation 
function has effectively a single channel. In particular, it can be applied to the calculation 
of special correlators involving only primary fields. Examples of such functions containing 
only Virasoro or parafermionic primary fields are presented below. 



2. Correlators of Virasoro degenerate primary fields 

We will first consider the following correlation function of Virasoro degenerate primary 
fields: 

(01 2 (^l) • ' • (pl2(z n )<Pl,n+l(Zn+l)} (2.1) 

Given the fusion rule [1] 

012 X (j) lr = 01, r _l + (pl, r +l (2.2) 

we see that the insertion of the field (pi, n +i in the last position effectively selects a single 
channel, i.e., a single term contributes from each OPE. This is most easily seen for the 
case n = 2 corresponding to the three-point function: substituting 4>i 2 x 4>i 2 = 0n + </>i3 
in the correlator and using the orthogonality condition 

{(f>lr{zi)(f>ls{z2)) = 5 r;S (2.3) 

it is clear that it is only the ^13 term does contribute. 
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Recall that rs , with r, s integers, refers to a Virasoro primary field that contains a 
singular descendant (hence the qualitative 'degenerate') at level rs and whose dimension 
h rs is linked to the central charge via a parameter t as follows: 



h rs (t) = \{r 2 - l)t + - A {s 2 -l)- t - l -{rs - 1) 
c(t) = 13-6 (t+ j 



(2.4) 



It is convenient to set 

« = i (2.5) 
In this notation, the dimension h\ n takes the form 

h ln = ±(n 2 -l)a-^(n-l) (2.6) 

so that the power of (z — w) of the leading term of the family <f> lm in the OPE of <f> lr x <j) ls 
is 

hir + h l8 - h lm = -(r 2 + s 2 - m 2 - l)a - -(r + s - m - 1) (2.7) 
In order to lighten further the notation, we will set 

01,r+l = V r (2.8) 

and define the structure constants c rs as follows 

V r (z)V s (w) ~ c rs V r+s + ■■■ (2.9) 

Therefore, in terms of the minimal-model structure constants, the c rs are 

c rs = C(i ir+ i) ) (i iS+ i) (1 ' r+s+1) = C(i )r+1 ) ) (i 5S+ i) 5 (i )T . +s+ i) (2.10) 

the last expression being symmetric with respect to the interchange of any two pairs of 
the three indices. 

Let us evaluate the four-point function 

(M*l)M*2)012(* 3 )014(*4)) = (Pl(*l)Pl(*2)Pl(*3)P 3 (*4)) (2-11) 

by the residue method. The following function 

F( Zl ) = z^z^z\l- x {Vx (^1)^1(^2)^1(^3)^3(^4)) (2.12) 
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turns out to an analytic function of z\ (there are no poles). Actually, it is simply a constant, 
as the behavior as zi^oo indicates, i.e., 

F( Zl ) ~ z^zr^r'^T (*i->oo) (2.13) 

where the last term is the contribution of the correlator per se, in which only V\{z\) 
contributes: V\{z\) ~ z^ 2hl2 . This constant can be evaluated in many different ways and 
in particular, in the limit Z\^>z 4 : 

lim F( Zl ) = lim z - 2 a z - 3 a 42- 1 ^ T (V 2 (z 2 )V 1 (z 3 )V 3 (z 4 )) 

1™ -v— a -a5a— 1 C H C 12 
= 2 ™/l2^13^14 5 a-l -a 4a-l 4a-l (2.14) 

i 4 z 14 z 23 z 24 z 34 

C11C12 



—a 5a— 1 5a— 1 

z 23 z 24 z 34 



Therefore, we have 



(V^V^V^Vsiz,)) = ._^_^_ a T!L 5a -i5 fl - 1 (2-15) 



z 12 z 13 Z 23 z 14 z 24 z 34 



This computation can be easily generalized to the case where there is an arbitrary 
number of V\ fields projected onto an appropriate V n field enforcing the single-channel 
constraint: 

(V 1 (z 1 )---V 1 (z n )V n (z n+1 ))=C 11 C 12 ---C 1 , n - 1 [J ^ II (n+2)a-l ( 2 ' 16 ) 

l<i<j<n Zi J l<£<n Z £,n+1 

It is not difficult to verify that this correlator is solution of the l2 singular-vector dif- 
ferential equation. Having computed our correlator without resorting to this differential 
equation, one could ask where does the singular nature of 0i2 enters, if it does at all. It is 
actually used right at the beginning, in specifying the fusion rules. 

Still using the residue method, we can derive the following other generalization of the 
four-point function (2.15): 

(Vl(zi)V n (z2)Vn'(z3)V n+n ' + l(z 4 )) 

_ Cl,n-\-n'Cn,n' 

a (3+n+n')a — 1 — nn ' a (n 2 +3n+nn')a—n (n' 2 +3n' +nn')a— n' 

(2.17) 



-na „— n 



Z \2 Z i 3 Z 1A Z 23 Z 2A Z M 



A different ordering in the evaluation of the constant representing the intermediate mero- 
morphic function forces the relation 



Cl,n+n'Cn,n' ~ Cl,nC n +l,n' (2-18) 

The solution of this recursion relation reads 

Cn,n> = C M"-<W-1 (219) 

Cl,l • • • Cx, n -l Cij • • • Ci jH '_x 

which, as it should, is symmetric in both indices. 

We now want to stress that the factorized expression (2.18) was actually coded in the 
correlator (2.16). This correlator was computed by contracting all the V\ fields from left 
to right. Equivalently, we could have stop this process at the m-th one and contract the 
remaining Pi's from right to left up to the m'-th. That yields (with n = m + m'): 

[cxic 12 • ••C l m _i][c 1 iCi2 ~'C\ 

(2.20) 

~ [cilCi2 • • • Ci >m _i] [C11C12 • • • Cl >m ' -l\C m ^ m i 

The comparison between the two results yields directly (2.19). 

Equation (2.19) shows that all the constants c n ^ n i can be calculated in terms of the 
ci >ri 's only.0 In order to calculate the structure constants ci jn , we need to evaluate the 
four-point function (ViV n ViV n ); but in this case, the two channels in the OPE V\ x V n 
do contribute, which invalidates the applicability of the residue method. This correlator 
can be evaluated by using the singular-vector equations [1] or by using screening operators 
[5, 6, 7]. But this will not be reconsidered here. Our main point was to unravel the 
factorization (2.19), as well as illustrating the residue method. 

We can similarly write down rather directly the expression for all correlators of the 

form 

(V ri [zx)V r2 (z 2 ) ■■■V rn (z n )V R (z n+1 )) 

TT 1 TT 1_ (2.21) 

C si,r 2 C S2 ,r 3 • ■ ■ Cs n -i,r n J_ J_ -r.rja H ri\(R+< 



l<i<j<n Z ij l<£<n Z £,n+1 



,[(R+2)a-l] 



1 We stress that the c„ in /'s form a particular class of structure constants, first in that they 
pertain to the restricted {4>i,r} algebra and second, because they are those in front of the 'highest 
field', i.e., the maximal value of m, in the OPE V r X V s ~ Pm- 
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with Si = J2 l j=i r j and Yh=i r i = R - 

Still more generally, we can also use the residue method to evaluate the correlators 

(021 (Zl) • • •02l(Zr-l)012(Zr) ' ' ' 012 (^r+s-l)0rs (^r+s)> (2-22) 

since they also involve a single channel. The corresponding structures constants satisfy 

(~i (r+r' — 1,8+s' — 1) _ /~i (r+r' — 1,1) (l,s+s — 1) _ „ /r> r>o\ 

^(r,s),(r',s') — ^(r,l),(r',l) ^(l,s),(l,s,) — C rr 'C ss ' 

Are these results completely surprising? From the Coulomb-gas representation point 
of view, these expressions are somewhat trivial: these are the very correlators that do not 
require the insertion of even a single screening operator (and for this reason they have 
not been considered in [5,6]). A Coulomb-gas correlation function without screening is 
simply that of a collection of vertex operators, the result of which being quite simple and 
well-known (see e.g., eq. (9.9) of [3]). All the Zij dependence of the correlators is recovered 
in this way, the remaining factors being simply the structure constants. 

However we stress that in our computation we do not require the free-field representa- 
tion. In this sense the present derivation is thus more fundamental, in spite of the fact that 
it is applicable to a rather limited class of correlators. It also reveals a simple factorization 
of some of the structure constants that may not have been obvious from other points of 
view, but whose generality is by now quite transparent. 



3. Correlators in parafermionic models 

3.1. Reviewing the Zjt parafermionic algebra 

The Zfc parafermionic conformal algebra [4,8] (see also [9,10] is generated by k 
conserved holomorphic (and similar anti-holomorphic) parafermionic fields ip n , n = 
0,1, - • • ,k — 1, with ifjQ = I and i/j^ = if^k-m with conformal dimension satisfying 
h"4> n = h-tpk-n- The set of conformal dimensions {hip n } is an input of a parafermionic 
theory and for the model, they read 

n(k-n) 

- J, 
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The ip n are primary fields that form a closed algebra specified by the following OPEs: 

n > (if) ~ - '\2nn' Jk ^ n + n ' 0^) (n -\- n < /c) 

Mz)^l(w) ~ (z _ ^^FwA V^+n-n'H (n + n' < k) (3.2) 



(2 _ w yn{k-n)/k 

where the central charge is fixed by associativity to be 

_ 2(fc-l) 
°~ (k + 2) 

The remaining OPE are obtained by conjugation and the condition 



(3.3) 



Cn,n' — Ck—n,k — n' (3-4) 

which implies, in particular, that c nj k-n' = c~k-n,n'- The constants c n ,ri are assumed to 
be real. They are fixed by the associativity conditions: 



2 r(n + n' + l)r(fc - n + l)r(fc - n' + 1) 



c n,n' 



T(n + l)r(n / + l)r(fc - n - ra' + l)r(fc + 1) 
Note that this yields c n ,fc-n = 1- They satisfy 



(n + ri < k) (3.5) 



Cn,n' — Cn,k—n—n' — Cn' ,k—n—n' (3-6) 

which reflects the symmetry of the three-point function: 

(lpn(zi)lpn'(z2)lpl +n/ (z3)) = 2nn , /k 2 n' -2n' (n+n'^k 2n-2n(n+n') / k ^.7) 

Z 12 Z 23 Z 13 

3.2. Zfc parafermionic correlators and their structure constants 

As a simple illustrative example of the application of the residue method to a 
parafermionic correlator, let us first rederive the well-known expression for the three-point 
correlation function 

<V'i(*i)V'i(*2)v4(*3)> (3.8) 

Given the structure of the OPEs, we know that the function 

F( Zl ) = £t k z£ /k (M*i)M22)4M) (3-9) 
11 



is a meromorphic function of z\ with a double pole at Z3 and no regular terms. Therefore, 
this function has to be of the form: 



F(zi) 



Cl2 d\ 



z 



13 



2l3 



(3.10) 



The first coefficient is 



a 2 = lim Zi 3 z% k 'z^ ,k '(V'l^OV'it^OV'fe)) 



(3.11) 



Since we will be interested also in the second order term, we will need to consider also the 
subleading contribution in the limiting value of the correlation function. Recall that in the 
OPE of (p a and <ph, the first two contributing terms in the conformal family of <p c are 



(j) a (z)(j) b {w) 



C a bc 



^ z — yj^h a +h b — h c 

In the present case, we need 

V>i(*i)V4(*3) ~ - 



<t>c{w) H — (z - w) d(p c (w) + ■ 



(3.12) 



Cl,l 
2-4/fc 
2 13 



4(^3) + j^jd4(z 3 ) + ■ ■ ■ 



(3.13) 



Therefore 



ci i(— l)~ 2 / h 

lim (^1(^1)^1(^2)^2(^3)) = ^im 2-4/fc 2-2/fc 



2:1^2:3 



2 z 13 
k z 23 



z 13 z 23 

Only the leading term contributes to the first coefficient, which is thus 

CIA 



(3.14) 



a 2 = 



2-4/fc 

2 23 



(3.15) 



The other coefficient is given by 



cii = lim 



d 



zi- 



z 3 dz\ 



2 2/fc -4/fc 
Z \3 Z \2 Z 13 



(3.16) 



It is simple to check that it is equal to zero. F{z\) has thus a single term; the expression 
of the correlation function under consideration is then 



(^1(^1)^1(^2)^2 (*&)) = 



This is the correct three-point function. 



c i,i 



2/fc 2-4/fc 2-4/fc 
Z 23 



(3.17) 



z 12 z 13 
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This computation can easily be generalized to the case where there is an arbitrary 
number of ijji factors : 

(Mzi)---Mzn)^ n (z n+ i))= n % n -i^jk ( 3 - 18 ) 

l<i<j<n Z ij l<e<n Z £,n+1 

Proceeding in a similar way, we can also compute the following correlator: 



Z 12 Z 13 Z 23 



1 

X 



(3.19) 



2-2(n+n' + l)/k 2n(l-(n+n' + l)/k) 2n' (1 - (n+n' + l)/fe) 
Z 14 Z 2A z 34 



We now use this last expression to extract a first result on the structure constants, by 
comparing the leading contribution of the limit zi—>Z2 of the correlator: 



(3.20) 



Z 12 



2n/k 2n'(n+l)/k J2(n+l)(l-{n+n'+l)/k) J2ri (l-(n+n' +1) /k) 



-12 ^23 ^24 ^34 

with the same limit calculated directly from (3.19). This gives: 

Cl,nCn+l,n' = Cn,n' Cl,n+n' (3-21) 

precisely the same relation found for the minimal models and thus whose solution is again 
exactly of the form (2.19). The observation of this factorization seems to have first been 
made in [11]. As already hinted at, it is typical of monomial (i.e., single term) correlators. 

In order to get the complete expression for the coefficient c n , n ', we simply need to 
evaluate ci i7l - This can also be done by the residue method.@ For this we need to consider 
the correlation function 

G' = (Mz^izzWUz^Uzi)} (3.22) 
with n > 1. As a function of z±, 

F'( Z1 ) = 4 n/ ^- 2/fc %4 2n/fc ^i^i)^(^)^(^3)^(^)) (3.23) 



2 That was not the case for Virasoro primary-field correlators since not all V\ correlators do 
have a single channel. The difference here is that the ip n 's have a single OPE channel. 
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is a meromorphic function of z\ with double poles at Z3 and Z4 and vanishing at infinity, 
so that 

w/ \ a,2 a\ b 2 hi 

F (*0 = 72- + — + 72- + — ( 3 - 24 ) 

z 13 z 13 z 14 z 13 

A straightforward computation yields: 

2n/k r, 

z 23 In Z24 

0/2 = 2n/k 2n-2n 2 /k «1 = ~0>2 ~~T~ 

z 34 z 24 Z 23Z34 



, __ vyt-i , , 2 ^23 

2 _ ^,- 2w /^2«/* 2(n-l)-2nV* X " 2 (fc - n + 1) ^34 



(3.25) 



-l,n-l 

,/fc 2(-, 
^23 ^"34 ^24 

from which we obtain: 

2n/k 



G' 



( Z14Z23 \ 1 | 1 _2n ^13^24 



V ^12^34/ zl 3 2/k zl™ 2n2/k I k Z23Z34 



, -2 Z 13 Z 24 , 2 ^14^23 "\ 

+ c l,n-l 9 „2 1 + 



1 ^23^14 V fe - Tl + 1 ^24^34 / 



(3.26) 



To extract the value of ci >n _i, we compare the expression of the correlator evaluated 
directly in the limit £12— >0 and 2 34 ^0, which yields 

hm 2 (M Z 1)M Z 2)4( Z S)^)) - 2n/k 2^-2^^ ( 3 - 27 ) 
z 3^ z 4 z 12 z 34 2 24 

with G" given in (3.26), evaluated in the same limit. One finds that the leading terms do 
not match, G' being more singular. Hence, the leading coefficient in G' must cancel, which 
forces: 

2 _ n(k-n + l) 

C l,n-1 - J, W-^) 

The substitution of this result into (2.19) leads to (3.5). Note, on the other hand, that the 
central charge is fixed by the subleading term of the correlator (ipi(z\)ipi(z2)ipl(z3)'ipl(z4)) . 



3.3. Some spin-field correlators and a relation between their structure constants 

A similar analysis can be applied to spin-field correlation functions. Here we use the 
notation cr^ for the holomorphic part of i-th spin field in the parafermionic theory, with 
o~o = crk = I. The fusion rules are 

01 x Oj = a i+j H (3.29) 



14 



where the dots stand for fields that are not parafermionic primary in the sense that their 
charge is not linked to their dimension as for the spin field, for which they are respectively 

hi = ~ £ \ and q e = I (3.30) 
2k(k + 2) H y ' 

These fusion rules are obtained from the coset realization su(2)k/u(l) [9]. Let 2jj be the 
structure constant in front of the term which is given explicitly in (3.29), namely 

(a n (z 1 )a n/ (z 2 )ai +n ,(z 3 )} = ^ n , bJbn{k _ n % ^ (fc - n _ n0 ( 3 - 31 ) 

where 



^12 ^13 ^23 



1 

(3.32) 



k(k + 2 



For the general correlator of n a\(z\) fields with erf, the residue method leads to the 
following simple expression: 

(o 1 (z 1 )---o 1 (z n )al(z n+1 )) = c 11 c 12 ---c hn - 1 Y[ -j- Yl (fc _ n)b (3-33) 

l<i<j<n H l<i<n Z i,n+1 



This again leads to the relation (2.19) - in tilde version - for the structure con- 
stants. Obviously, this relation could be derived as previously, by considering the correla- 
tor (cricr n cr n /cr^ , n , , j), which leads to (2.18) whose solution is the tilde version of (2.19). 
It is simple to check that the explicit expression found in [4] for the c n ,ra / S 

~ 2 = liX)l{n + n' + l) 7 (fc - n + l) 7 (fc - n' + 1) 

Cn ' n ' 7 (n + l) 7 (n' + l) 7 (/c - n - n' + l) 7 (Jfe + 1) 1 J 

with 

7(«) = r (^) (3.35) 

indeed satisfies the factorization (2.19).i This is at the roots of the curious similarity noted 
in [4] between the parafermionic structure constants and the spin-field ones. 



3 Note that the factorization by itself cannot lead us to the above expression: the exact value 
of ci, n is required and it can be obtained through the calculation of the correlator (a\o- n o'lcT} l ) . 
But this cannot be evaluated by means of the residue method because there are two contribut- 
ing channels. In [4] it was evaluated through the coset representation. The coset approach is 
usually not a convenient way of computing correlation functions - (cf. [3], chap. 18); however, 
the simplicity of the present coset, namely su(2)k/u(l), allows for a direct factorization of the 
correlators into a WZW piece and a free boson one. The correlator could also be calculated from 
the parafermionic singular conditions. This calculation will be reported elsewhere. 
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3.4- A Zl ' correlator 



As already pointed out, in a parafermionic model, the conformal dimension of the 
parafermionic fields is an input. It is however constrained by the Z& invariance itself, which 
requires h 1 p n = h^t and the monodromy invariance of the correlators. These constraints 
are satisfied for 

= (3.36) 

for any positive integer f3. However, the underlying associativity conditions have to be 
checked anew for each value of f3. The defining OPEs are given by (3.2) but with (z — w) 
replaced by (z — w) 13 and c n ,ri — *<^ n " ^ n ^h e following, we denote the corresponding 
parafermionic model as Z^ , with Z& = Z^ 1 "*. The results of the associativity conditions 

(2) 

for the case (3 = 2 are been presented in [4] (cf. their appendix A) and the special Z 3 
model is studied in [12]. (We will report elsewhere on a detailed analysis of the Zj^ 
models.) 

Let us consider a sample correlator of the Z^ model, namely (ipiipiiplipl) , from which 
we construct the meromorphic function 

F( Zl ) = zti k z; 3 4/k z; 4 4/k (Mzi)Mz2)4(^l(^)) (3.37) 

It has poles or order 4 at 23 and 24 and behaves as l/zf as 21^00, meaning that there is 
no analytic piece. The computation of the two principal parts requires thus the knowledge 
of the first three subleading terms in the OPE ipi(z)ip\(w) . That makes the computation 
much more involved than in the (3 = 1 case. However, there is a simple trick that allows 
us to avoid going so deeply inside the conformal block. Since F ~ 1/zf, we can reduce 
the order of the two poles by two, at the price of adding a constant term, by multiplying 
F(z\) by z\ 3 z\ 4 . The transformed meromorphic function reads thus 

F'( Zl ) = z 2 13 zlF( Zl ) = ^ + ^± + ^ + ^ + Co (3.38) 

z 13 z 13 Z 14 Z14 

For this computation we require the knowledge of only the first subleading term in the 
OPE ipi(z)ip\(w). A direct but rather long calculation gives the following expression for 
the correlator: 

2 



(2)V 
11 ) 



+ 



2 34 4 Z 12 Z 34 



(z 12 Z34) 4/k (z 1 sZ 14: Z23Z24) 2 4 / k \ \ ) ^23^24 k 2l42l 3 2 2 3^24 

. ^34^24 / \ ^34^23 / -> 

+ — 2~\ Zl2Zl4 ~ ^13^24 J H 2 2~~ 1^14^23 — ^12^13j 

Z 14 Z 23 Z 24 Z 13 
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(3.39) 



If we compare this expression in the limit Zi^zs with the correlator evaluated directly in 
this same limit, we find that for c expressed in terms of a parameter A as 

4(fc-l)A(fc + A-l) 

(£; + 2A)£; + 2A-2) 1 ' 

(2) 

the structure constant c-q becomes 

({2) y 2(fc-l)(A + l)(fc + A-2) 



k\(k + X-l) 

in agreement with [12]. To check that the 7L\ ' central charge is unconstrained, we need to 
compare the sub leading terms. 

(2) 

Here again the constants c n n , have the factorization property (2.19), so that all con- 

(2) 

stants can be expressed in terms of c\ n , whose determination requires the evaluation of 
(V'lV'nV'iV'n)- Again the calculation can be reduced to the evaluation of poles of order two. 

(2) (2) 

The expression of all the structure constants c n n , are given in [12]. Curiously, c n n , has 
another factorization, namely as 

= Cn,n'C?n,n'(A) , with d n , n >{l) = C n , n ' = C^, (3.42) 

whose origin appears somewhat mysterious. 

The structure constants for the models Z ( jf >2 ' > can be computed in the same manner. 
Yet there exits no results concerning these theories. 

Coming back to methodological aspects, with the computation presented in this sub- 
section, we wanted to emphasis that the evaluation of the principal part can be sub- 
stantially simplified by an appropriate modification of the integer powers of the prefactor 
multiplying the correlator. Here, poles of order 4 have been transformed into poles of order 
2. This is not a purely technical issue. The determination of higher sub leading terms in 
the OPE ip^ requires the knowledge of the chiral algebra underlying the Zj^ models. More 
precisely, the required information can be extracted from various associativity constraints 
but it is clear that knowing at least the gross features of the underlying chiral algebra is 
useful. For the Z^ model, it is the WAk-i algebra. However, the chiral algebra of the 

(2) 

Zj. ' models is not known. 
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4. Conclusion 

We have thus exposed, in its full generality, a method for calculating special CFT 
correlators based an a direct application of the residue calculus. Let us first summarize 
the method and restate the limits of its applicability. 

We consider the correlator (A(z)Y\iBi(zi)) as a function of z. The first step is to 
transform it into a meromorphic function of z. This is already guaranteed if the OPEs 
A(z)Bi(zi) are Laurent series, which requires that the conformal dimensions of all the fields 
that appear in the singular terms of the OPE A(z)Bi(zi), as well as A and Bi, differ by 
integers. This situation pertains to correlators of chiral-algebra generators (in which case 
all the fields have integer dimension). However, when A is not a symmetry generator, the 
above condition is rarely verified. Generically, the OPE A(z)Bi(zi) contains a number of 
channels associated to primary fields whose dimensions do not differ by integers. In such a 
case, one can still construct a meromorphic function if either there is a single channel, or 
more generally, if a single channel, say C i: contributes to the correlation function. This, 
however, does not ensure that Yla + — is integer, i.e., that the OPE projected onto 
the Ci channel is a Laurent series. But the cure at this point is simple: one multiplies the 
OPE by a fractional power of z — Zi suitably chosen to eliminate the algebraic singularity, 
leaving thus pole-type singular terms. In other words, when the effective single-channel 
requirement is satisfied, the correlator can be transformed into a meromorphic function of 
z by multiplying it by an appropriate factor Yli( z ~ z i) di • This meromorphic function is 
then expanded in partial fractions, that is, as the sum of the principal parts at the different 
Zi, including possibly the principal part at infinity. 

This method is certainly not new and we try in appendix B to trace it back in the 
literature. However, its explicit spelling out as well as the identification of its inherent lim- 
itations appear to be new. The method has been illustrated here with various examples, 
including some correlators that do not involve conserved currents. As a practical applica- 
tion, we have worked out a detailed derivation of the parafermionic structure constants. 

Appendix A. A three-point function computed from the infinite series 

In this appendix, we consider the calculation of the three-point function {T{z\)T{z2)T{zz)) 
using the brute force infinite-series method, where one OPE is replaced by its full infinite 
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series and show how this series can be summed exactly to reproduce the simple three-point 
function expression. 

The long road computation will be compared to the few step computations that results 
from the application of the conformal Ward identity: 



(T( Zl )T(z 2 )T(z 3 )) = (T( Zl )T(z 2 )T(zs)) + (T(z 2 ) T( Zl )T(z 3 )) 

eft + 2T^ + a^N V + / + 2^) + 

Z 12 Z 12 Z 12 J I \ V^13 z 13 *13 

c 2c c 2c c 

z \l z 2Z Z 12 Z 23 z 13 z 23 Z 13% 23 Z 1 2 Z 1 3 Z 23 



(A.l) 



Let us now turn to the expression that results from the replacement of the OPE 
T(zi)T(z 2 ) by its infinite series: 



T( Zl )T(z 2 ) = r( Zl )T(z 2 ) + (T( Zl )T(z 2 )) 

\*12 Z 12 Z 12 J ^ n! 

where T^ n ) = <9 n T and {T{z\)T{z 2 )) stands for the normal ordering of the product 
T( Zl )T(z 2 ), defined as 

{AB){z)= ^f^z A{x)B{z) (A - 3) 
The substitution of (A. 2) into the three-point function yields 



2c A zk 



(T( Zl )T(z 2 )T(z 3 )) = r + J2 ^((TT^)( Zl )T(z 3 )) (A.4) 

z 12 z 23 z 12 z% 3 ^ n\ 

In the last three-point function, we have to find the term proportional to c Zi 3 ~ n which is 
the only contributing one; it is obtained by standard methods (see e.g., [3]): in 

T(z 3 ) (TT^)( Zl ) = ^-j | (rte)T(x)j T^( Zl )+T(x) (t(z 3 )T^( Zi ) 

(A.5) 

only the first piece contributes to cz 3 ±~ n and its different contributions add up (6 + n)(3 + 
n)!/12. Therefore, we have 

ifrpf xrpf \ \ rri / u c ^ (3 + n)!(6 + n) z? 2 
({T{zi)T(z 2 ))T(z 3 )) = — ~\ ( A - 6 ) 



n=0 '"• Z 13 
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The infinite series can be summed as follows: 





(A.7) 



+ 



24zf 3 z 12 



23 



so that 



((T( Z1 )T(Z 2 ))T(Z S )) = -2^(3^23 + 2Z 12 ) = ~^(Z23 + 2z 13 ) 



(A.8) 




The substitution of this expression into (A. 4) reproduces (A.l). 

Appendix B. Correlators as meromorphic functions through the literature 

As applied to correlators containing symmetry generators, the residue method could 
be traced back to the pioneer work of Belavin, Polyakov and Zamolodchikov [1]. These 
authors explicitly consider the correlator (T(z) f|j 4>i{ z i)) as a meromorphic function of z 
with poles at Zi whose residues are fixed by the conformal properties of the field (pi(zi) 
(cf. the discussion after their eqs (2.8) and (3.3)). The method is also used in [2]. In 
this seminal paper, Zamoldchikov has launched the exploration of extended conformal 
algebras through the study of the associativity conditions of a number of cases containing 
a single extra symmetry generator. Correlation functions are computed by considering 
only the singular terms in the OPEs. There is again an explicit reference to the residue 
calculus (cf. eqs (2.1)-(2.4)) that suggests an underlying complex-analysis interpretation 
of the exposed computations. But note that such an approach is not mandatory since the 
calculation method could be justified by means of the Ward identities associated to these 
extra conserved currents. 

Meromorphicity is the central theme of Goddard's proposed formalization of conformal 
field theory [13]. However, the meromorphic point of view here is not implemented at 
the level of computing correlation functions. In this context, meromorphicity is used to 
establish locality which in turn becomes the corner-stone property for the construction of 
the conformal field theory. This formal procedure has been much developed in a sequel 
work with Gaberdiel [14]. Here the emphasis is placed on the reconstruction the whole 
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theory out of meromrophic amplitudes, via an approach inspired by the early works in 
dual models where the space of states was originally built out of the conjectured dual 
amplitudes. In this very paper however, the residue method is explicitly invoked for the 
calculation of some correlators involving symmetry generators. In that regard, they present 
some results of Frenkel and Zhu [15], which have devised a nice combinatorial method 
for handling such correlators (involving an arbitrary number of either T factors or affine 
Lie algebra generators - summarized in sections 5(b) and 5(c) of [14]). But again, the 
residue method is not systematized and its applications are restricted to the elimination 
of conserved currents in correlation functions. 

It is in the context of a non-meromorphic theory that the partial fraction expansion 
technique has been mentioned in the most explicit way, namely as a natural tool for eval- 
uating the parafermionic correlator (ipx ■ ■ • V'lV'i ' ' 'V4) [4] (cf. the discussion between eqs 
(3.9)-(3.12)). In this paper, the authors also give the value of the parafermionic structure 
constants. Even though there are no indication concerning the way these have been com- 
puted, it is natural to guess that their calculations have been done roughly along the lines 
presented here. 

A variant of the partial fraction expansion method has been used explicitly in [11] to 
work out in detail the associativity conditions of the Z^ 2 ) models (cf. their section 3B). The 
approach used there is superficially a little more complicated that the one presented here 
in that the four-point correlation functions are transformed into meromorphic functions 
(actually, into polynomials) of the cross ratio instead of functions of the position of one 
field. This procedure prevents an immediate generalization to higher-point functions. 

As applied to correlators that do not contain symmetry generators, we found after- 
wards a single reference to the residue method: this is in appendix E of [1] (cf. eqs 
(E.14)-(E.17)). It is used there to calculate the Ising correlator (ip(z)a(zi) ■ ■ ■ h{z2m)) 
where tfj is the Ising fermion and a and \i are respectively the spin and disorder fields. The 
correlator is transformed into a meromrophic function of z as follows: 




(B.l) 



Actually, all singularities in z are removed by this transformation since 



tp(z)a(w) 



a(w) 



(B.2) 



(z-w) 1 / 2 



(z-w) 1 / 2 
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Observe that these OPEs have a single channel. F(z) is thus an analytic function. The 
form of this analytic function is determined by its behavior at infinity. Since F(z) ~ z M ~ x 
as z^oo, it is necessarily a polynomial of order M — 1, which can be written as 

M-l 

F(z) = J2( z - Z ^) k 9k{zi) (B.3) 
fe=i 

(disregarding the dependence upon the anti-holomorphic variables) .@ The coefficients gk 
are determined by enforcing the correlator to be a solution of the singular- vector differential 
equation of the free fermion (= </> 2 i) . Quite interestingly, the coefficient go gives the value 
of the correlation function without the fermion. 

More recently, Dotsenko [16] has devised a nice way of handling the parafermionic 
computations by a method which is close in spirit to the computation of [1] just described, 
hence to the partial fraction expansion. Applied to four-point functions with three points 
fixed at the special values 0, 1, oo, the idea is to factor out the branching or poles singu- 
larities of the resulting z function and then determine the remaining polynomial in z that 
completes the correlator by considering successively the correlator in the limits where z 
approaches the three fixed points. A detailed application of this method is presented in 
the following appendix. 



Appendix C. The Dotsenko's method to test associativity 

We will illustrate the method initiated by Dotsenko by reconsidering the 
parafermionic correlator (tpi(zi)ip n (z2)ipl(zs)tpl l (z4)) . We fix three points at the standard 
values 0, 1 and oo and use the following convention for a 'prime correlator': 

(XA(oo))' = lim z 2 n hA (XA(z n )) (C.l) 

In particular, we have 

(M^Uoo))' = i , (Mz)M^ n+n/ (oo)Y = {z _°^ nnl/k (C2) 

The choice of the expansion variable, here z — Z2M, is conventional; another Zj could have be 
chosen instead of zim and the power series could also have been written in powers of z simply. 
Writing the expansion in terms of a difference between two variables has the advantage of taking 
care of translation invariance. Note that there is a misprint on the value of the upper limit of the 
sum in (E.16) of [1]. 
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We consider thus (ipi (0)ip n (z)ip\ (l)^(oo)). Regarded as a function of z, this has singu- 
larities at and 1, that is 

G{z) = (Mo)M*)4Mti{°°)Y = z2n/k{ ^%- 2n/k (C3) 

(we do not care about the phases which all cancel at the end) . P n is a polynomial of order 
n. In other words, by multiplying the correlator by the prefactor z 2n l k {z — i) 2_2n / fc ; we 
transform it into a meromorphic function, which turns out to be analytic. 

The order of P n is fixed by considering the limit where z-^oo, where we readily see 
that (/(oo) = 1, i.e., 



lim Q ~ lim (M0)4W)(M^n(oo)Y = 1 (C.4) 
This implies that n = 2 and, in addition, that the coefficient of the term z 2 is 1, that is 

P 2 (z) =a + a 1 z + z 2 (C.5) 

Thus, in order to completely determine the correlator, we only have to fix these two 
constants. 

Consider first the limit where z— >0, keeping track of the first subleading term. The 
correlator Q(z) becomes then 

C tn ^ 2(fe-n-l)z | 

z 2n/fc( z _ 1 )2-2(n+l)/fc \ ^ fc(n + l)(z - 1) J 

^2 



ci.n f 1+ 2fc(fc-n-l)z 



^2n/fc |^ /.( n + J) 

This is to be compared with the expansion of the rhs of (C.3): 



lim 



a + a lZ + z 2 1 1 / 2(k-n) 

a + z ai + a - 



2 ^ [ Z 2n / k (z - l)2-2n/fc J ^2n/fc 

This yields 



fc 



(C.7) 



-2 2a / n c n 

a = c 1>n , a i = -^ry ( c - 8 ) 
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Next we consider the limit z— >1. Here it will suffices to keep only the leading term. 
We have thus 

which is to be compared to 

f a + ai^ + z 2 ] ap + ai + l , , 

z ™ | ^n/fc^ _ 1 )2-2n/fc J ~ (z - 1)2-2n/fc 

That forces 

a + ai + 1 = c\^_ x (C.10) 
which together with (C.8) yields the recursion relation 

^-i='1. n (^)+l (an) 

This relation can be solved without knowing the explicit value of en but simply by enforcing 
the condition 

Cl, n = Ci,fc-n-l (C12) 

This readily implies 

c\ k _ n = c\ k _ n _ x + 1 (C13) 

On the other hand, replacing n by — n in (C.ll) yields 



Eliminating c\ k _ n from the last two equations lead to an algebraic relation for c\ k _ n _ 1 
which gives 

^ _=2 _ (n + l)(k-n) 

c l,fe-n-l — c l,n — r ^.10J 



and we recover the expression of ci jn obtained previously in (3.28). 

At first sight, it seems that the applicability of this method depends critically upon 
the fact that the field evaluated at z is the conjugate of the one at infinity. In that case, the 
other two fields are conjugate of each other, which ensures that the fractional powers of z 
cancel in the denominator as z^oo. That certainly ensures the polynomial character of P n 
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as in the large z behavior ip n (z) is simply projected onto 1/^(00) and {tp n (z)^(oo))' = 1. 
But consider instead the correlator 

(^(0)^(^)^(1)^(00))' = z2n/k{z Q ^ 2n _ 2nyk (C16) 

It is not clear at once that Q(z) is polynomial here. The point however is that the large z 
limit of (A(0)B(z)C(l)D(oo)}' is actually given by z hB+hl3 ~ hB / z 2tlB where E is the single 
contributing field appearing in the OPE of B and D. Indeed, by considering at first the 
limit 24^00, we wash out the contribution Z2f +h °~ hE that needs to be reinserted at this 
point in order to get the right large z behavior. The term z~ 2Hb simply corresponds to the 
large z behavior of the B field. Returning to our problem, we need to compare the large 
z limit of the rhs of (C.16) with z 2 ~ 2n+2n2 /k-2n/k ^ z 2 s hows that Q(z) is indeed 

a polynomial of degree 2. 

It should be clear from this example that the Dotsenko's method has the same intrinsic 
limitations as the partial fraction expansion described in the main part of the article. In 
particular, it generically applies to correlators that have a single contributing channel. 
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